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We present a new approximation method for solving the equations of motion for cosmological 
tensor perturbations in a Randall-Sundrum brane- world model of the type with one brane in a five- 
dimensional anti-de Sitter spacetime. This method avoids the problem of coordinate singularities 
inherent in some methods. At leading order, the zero- mode solution replicates the evolution of 
perturbations in a four-dimensional Friedmann-Robertson-Walker universe in the absence of any 

' tensor component to the matter perturbation on the brane. At next order, there is a mode-mixing 

' effect, although, importantly, the zero-mode does not source any other modes. 

o ' 
(N 

< 

' The concept of brane-worlds is not a very new one [l|, but has risen in popularity in the last five years. Central to 
' , this high level of interest has been the observation that the experimental bounds on the inverse square behaviour of the 
^ ' gravitational force are surprisingly large, with the current bound being about 10~'*m 0; by contrast, the equivalent 
I bound on the Coulomb law is about 10~^^m. In a seminal paper, Arkani-Hamed, Dimopoulos and Dvali Q observed 
I that, if the Standard Model interactions were confined to a four-dimensional subspace, the length-scale of extra 
' dimensions could be made as large as these bounds on Newton's law. A further revelation came from the work of 
^-H Randall and Sundrum 0, , who demonstrated that the extra dimensions could be made infinite in extent if one 
■ removed the assumption of a product geometry. 

, We will focus on the single-brane Randall-Sundrum (RS) model The cosmological background solution for this 
model was found shortly after the model was suggested ,6j; however, the question of how perturbations evolve in such 
a background is still unresolved. Some progress has been made, for example, it is possible to analyse the evolution of 
perturbations perturbations during a primordial de Sitter phase |3j 13 of the brane- world assuming an initial vacuum 
state. Our interest here is in the a more general cosmological background dominated by either radiation or matter. 

Several approaches have been proposed to study this problem. One is to try to construct an effective theory in four 
dimensions and study perturbations in that. This is particularly useful for brane- world models with a finitely sized 
extra dimension since it is then not possible to excite Kaluza-Klein states below a certain energy; however, it is less 
well motivated in the case of the single-brane model. In this method one ignores most of the higher-dimensional effects 
by reducing the problem to a scalar-tensor theory, so even where it is an accurate approximation it can never offer 
5J] ' a "smoking gun" for brane-world models by predicting a uniquely higher-dimensional effect. Other approaches have 
[ been to try to solve the full five-dimensional problem, either numerically or analytically, using some approximation 
scheme. Typically, these methods fall into two categories: brane-based and bulk- The brane-based methods use a 
Gaussian normal (GN) coordinate system, constructed from a congruence of geodesies normal to the brane, where the 
brane is chosen to lie at a fixed value of one of the coordinates. These GN coordinates make the boundary conditions 
at the brane simple at the expense of rendering the equation of motion for the perturbations complicated. The GN 
coordinates also suffer from a coordinate horizon, meaning that they do not cover the whole spacetime. In the case 
of a brane with the induced metric of de Sitter space, the coordinate singularity coincides with the horizon, and, as 
mentioned above, it is possible to solve the problem 

In the bulk-based approach one typically chooses a coordinate system that makes manifest the symmetries of the 
bulk, resulting in a simple equation of motion for the perturbations; the penalty for this is that the brane is not at 
a fixed value of one of the coordinates and, consequently, applying the boundary condition is complicated. In the 
special case of a brane with a Minkowski induced metric, the brane is at a fixed value of one of the coordinates and the 
perturbation equations can be solved; indeed, this was done by Randall and Sundrum did in their original paper 
One major advantage of the bulk-based approach is that it does not have any coordinate singularities. The relation 
between the GN and bulk coordinates is presented in Appendix 1X1 

Several methods for solving this problem have been suggested in the literature, all of which have some limitations. 
One possibility is to make an educated guess about the behaviour of the bulk effects 10] and treat these effects as 
a source term in the usual equations for the evolution of perturbations in four dimensions. The obvious drawback is 
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that the guess might not be educated enough; furthermore, the perturbations on the brane can source the bulk effects, 
meaning that a simple ansatz is unlikely to give a realistic answer. To go beyond this, many authors have suggested 
methods^_of approximating the problem, for example, the near-brane approximation [Tll | or the gradient expansion 
method 12]. The main limitation with both of these approaches is that they rely on the GN coordinates, which have 
a coordinate horizon. 

In this article, we present a new approximate method for solving the evolution equations for tensor perturbations. 
There are two stages to this: a change of coordinates to a hybrid system — a bulk-based system modified to make 
the boundary condition tractable within a well-motivated approximation scheme — and an adaptation of the method 
of multiple scales |l3l Il4| to find approximate solutions to the equations of motion. Intuitively, this approximation 
works by assuming that the boundary condition evolves slowly compared to the time-scale of the solutions. We will 
comment at the end on the consequences of the solutions we obtain and on how we hope to extend this method to 
work for scalar perturbations. 



II. THE EQUATIONS OF MOTION 

We will start from the Poincare coordinates, as advocated in Ref. |T5|. where the metric is manifestly conformally 
related to Minkowski space. The line clement has the form 

ds2 = ^T^^p dx"dx^ = ^ ( - dT^ + % dx'dx^ + dZ^) , (1) 

where / is the AdS length-scale. We have assumed spatial flatness of the background cosmology here. This coordinate 
system makes the five-dimensional linearized Einstein equations simple, the price being that the brane has locus given 
by Z = l/a, rather than being at a fixed value of one of the coordinates, making the boundary condition much more 
difficult to impose. If we write the perturbed line element as 

ds^ = ^ (vap + Kp) dx'^dx'^ , (2) 

and choose the transverse-traceless (TT) gauge, following Ref. where hap satisfies h^z — 0, rj^^haf) = and 
^ — 9''-i9°'tJ.9"i3^ph^v = 0, with g and V representing, respectively, the projected metric and derivative on the 
brane. This setup has the advantage that the resulting perturbed Einstein equations exactly soluble. After 
making a Fourier transform in the spatial directions parallel to the brane, the equation is simply 

d^h d^h 3 9/1,2, 

ar^-az^ + zaz + ^^-o, (3) 

where the indices on h have been dropped for convenience. This equation can be solved by separation of variables to 
give mode functions h„i = exp ( ± Wk^ + m^T) Z^B2 {mZ), where B2 is a Bessel function of order 2. There will be 
five different polarizations of the graviton in five-dimensions. However, in this gauge the position of the brane is not 
on the same locus as for the background but is displaced, an effect which has been dubbed "brane-bending" in the 
literature. This can be written as An^, that is, a displacement of A from the background brane position along the 
direction of the unit normal n, which is given by |15| 

n^ = i(£,0,0,0,-v/rT7^) , e^-L^^ — ^lH, (4) 
a \ / a dT a dt 

where a is the scale factor on the brane, and r and t are, respectively, conformal and proper time on the brane. 
Combined with the five polarizations of the gravitational perturbations in the bulk, this give six degrees of freedom, 
corresponding to the scalar, vector and tensor degrees of freedom in four-dimensional cosmology. Since the pertur- 
bation of the position of the brane is a scalar degree of freedom, we can ignore this effect when studying tensor 
perturbations to linear order. This approach should also apply to the bulk perturbations which give rise to vector 
modes on the brane, although we will not discuss these here. 

The boundary condition at the brane is given by the Israel junction condition >1^ . relating the jump in the normal 
derivative of the metric to the matter supported on the brane. (With reflection symmetry imposed between the two 
sides of the brane, the jump is equal to twice the value on one side.) For tensor perturbations this is simply 
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where the r.h.s. is the tensor part of the perturbation of the brane stress-energy-momentum tensor. For most of the 
history of the universe, there is no matter source for the tensor modes so the r.h.s. of (jSJ will reduce to zero. Only 
when the CMB photons develop a quadrupole moment in the late universe is this assumption no longer valid. 
In summary, the problem we must solve for the evolution of the tensor perturbations is 



Qf2 



a^/i 3 a/i 

'dZ^^'z'dZ^ ' subject to 



r dh 



Z=l/a 



dh 
df 



(6) 



Z=l/a 



which, although seemingly innocuous, is extremely difficult to do exactly. We will try to find approximate solutions 
based on the observation that there are two independent dimensionless parameters in the problem: e = IH and p = kl, 
which are, respectively, the Hubble parameter and the frequency of the perturbation measured in terms of the AdS 
length-scale. The presence of two scales means that the solutions should have time dependence on two different scales 
— much like an amplitude- or frequency- modulated radio broadcast, which is a signal of a hi gh frequency with its 
amplitude or frequency changing on a much longer time-scale. The method of multiple scales [ij, [Tjl is specifically 
designed to find approximate solutions in precisely these situations. 



III. NEW COORDINATE SYSTEM 



We will be constructing an approximation scheme valid when e = IH is small. If one naively applies the method 
of multiple scales as described in section IIVI one obtains an incorrect answer. This is because, to get the correct 
evolution at zeroth order, it is necessary for the boundary condition to have no first-order term. In order to make the 
method work, we devise a new set of coordinates designed to remove the first-order term in the boundary condition. 

From the change of coordinates in the appendix, we have that, on the brane, 

dr 



— = v/TT7^, (7) 

where r is conformal time for an observer on the brane moving with the Hubble flow. We can then evaluate 

da dr dt da a? r 

dr-dTd7d^--^^^- 

If we also assume that the Friedmann law — 8ttGp/3 holds to sufficient accuracy, that is, we are late enough in 
the history of the universe for the correction to be ignorable, then we can compute that 

_ = _4^G(p-fp) = --^i^^ (9) 

for a cosmological fluid with polytropic index F satisfying p + p = Tp. This gives 

d£ ar 



If we make the change of coordinates 



^=-^ae^yiTi^. (10) 



the equation of motion becomes 



d'^h d^h 3dh „ ( . d^h dh\ 



where p = kl is the dimensionless frequency of the mode. The boundary condition is simply 

0{e'). (13) 



dh 
5C 



= l/a 



Note that the brane is still at ^ 1/a in this new coordinate system. 
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IV. MULTIPLE-SCALES ANALYSIS 

For the purpose of multiple-scales analysis [13. H^ . we introduce two new time variables: "fast" time /, and "slow" 
time s, given by 

f = V, = (14) 

a 

Note that these are dimensionless (as is ^). The variables / and s are both, of course, functions of rj so the definition 
is a technical construct to allow us simply to handle the fact that the boundary condition changes with time on a 
scale related to s. Let us list our various time coordinates for clarity: 

T is the time coordinate in the static coordinates, 

7] is the time coordinate in the new coordinates, 

T is conformal time on the brane, 

t is proper time on the brane, 

/ is the "fast" time coordinate, 

s is the "slow" time coordinate. 

We assume that e is small, which is valid except in the very early universe. In the usual application of the method 
of multiple scale 0, the small parameter e is a constant; here it is a function of T. We then write derivatives 
w.r.t. 1] in terms of the fast and slow time derivatives as 

# = ij-^^+Oie'), (15) 



2 - 2£^+0(e^). (16) 



drj df ds 

drf 9/2 dfds 
Expressed in terms of these fast and slow variables, the equation of motion is 

subject to the boundary condition H13|l at ^ = s. 

V. SOLUTIONS 

We will use the separation of variables technique, where one tries to find a solution of the form 

h = ^ *„,(/, s,£)$,„(^,s,e) dm. (18) 

We will write ^'to(/, s, e) and 3'm(C, s, e) as series in e: 

*„,(/, s, £)« VI, W(/,s)+e*« (77, s) + ... , <J>™(f,s,£)«$,[^)(f,s)+£ci>W(^,s) + ... , (19) 

which we will require to be asymptotic as e — > 0, and we will require that the $m(^, s) satisfy the boundary condition 
(|13|l individually. 

A. Zeroth-order solution 

We will construct the solution order by order in e. The equation of motion to zeroth order is 

1 , 1 92$(o) 3 1 a$ 



^(0) Qp ^ $(0) Q^2 ^^(0) 

which is separable with solutions 

*^He,s) = |r<2(ri(ms) J2(mC)- Ji(ms)y2(m0) , (22) 
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with being chosen to satisfy the boundary condition H13|) to zeroth order: 



a$(0) 



= 0. 



(23) 



The zero-mode solution s) = s is the hniit of (|22|l as to ^ 0. Of course, one has the freedom to redistribute 

the s dependence between $ and \1/ so as to leave the product unchanged. We can use the relation derived from the 
Wronksian of Bessel functions, given in the Appendix B, to evaluate $o°^(Ci s) on the brane as 



1 



air) 



(24) 



Note that what would be constants, fixed by the initial conditions, in a standard solution of (|20() arc now functions 
of the slow time variable, s, and it is these functions which encode the crucial evolution of the amplitude. 



B. First-order solution 



The function Am{s) in expression (|21|l is undetermined by the zeroth-order calculation; in order to calculated it, 
we must consider the first order equation and construct the so-called secularity condition which forces the solution to 
asymptotic. The order e terms in the equation of motion (|17|l give us 




3 9$ 



(1) 



$(0) 



dz^ z dz 



+ to2$(i) 



dfds df ds s df 

and the boundary condition to order e is the same as at zeroth order 



$(0) + i 



0. 



We still have some freedom in constructing the solution, which we shall use to impose 



3(9$ 



(1) 



so that satisfies the same equation as ^m' . Given that the boundary condition for ^m' is also the same as $m , 
it is actually be possible to impose the stronger condition 
the orthogonality relation. 



.(1) 



(25) 



(26) 



(27) 



.(0) 



without any loss of generality. We can then use 
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$W $W = ^m(^(j^ims))^ + (yi(TOs))') S{m - n) , 



(28) 



to write the equation of motion as an equation for ^m-* , giving 



.2,7,(1) 



dp 



-^2 ^(1) 



2,T,(0) 



1 d-^ 



(0) 



dfds s df 



Ji(tos)2 + Yi(tos)' 



dn 



(0) 



df 



de<j,(o)5'i> 



(0) 



e " ds 



dn 



d-^ 



(0) 



df 



^2 »» 



(0) 



d^ 



(29) 



The calculation of the various integrals involving Bessel functions is explained in Appendix B. We can write the 
equation of motion for ^'m'' as 



.2,1,(1) 



~dp 



^(1) 



2,T,(0) 



92^ 

dfds 



d_ 
d^ 



logy' (Ji(ms))^ + (Yi(tos)) 



(0) 



df 



(0) 



dn A4j 



dm 
df ' 



(30) 



6 



for m 7^ 0, and 



2,T,(") 







dfds 



(0) 



for the zero-mode. The trace-free matrix Mmn is defined as 



Mr 



dn Mq, 



' df ' 



(0) 



(31) 



(32) 



for m ^ n and is calculated in Appendix B; these terms do not contribute to the secularity condition. 
Substituting in the solution (|21|) for \ we get terms on the r.h.s. of (|31|) for the form 

^ ds 



(33) 



This is a forcing term in equation H30|l which, as a function of the fast time /, has the same frequency as the natural 
frequency of the l.h.s. and so will cause a resonance. If £ decays as I/77 or slower, this is forbidden since the h'^^^ term 
would not be of asymptotic order. So we can impose the secularity condition that Aq{s) must satisfy 



ds 



(34) 



where the are constants independent of s. 

For the massive modes (m ^ 0), the requirement that (|30|l be free of forcing terms gives the secularity condition 



_d 
d^ 



-1/2 



= 



(35) 



where, again, the are constants independent of s. We will choose cx m so that the limit as m —s- is a 
constant. With the secularity condition imposed, equation of motion for ^'m^' is then 



dp 



which will have solution 



o o •' ' '■mn 



(36) 



(37) 



In order to determine Bm{s), one would have to consider its secularity condition, coming from the second-order 
equations. 



VI. MODE-MIXING AND THE ZERO-MODE 



A great advantage of the method of multiple scales is that the leading-order solution includes the slow-time variation 
derived from the secularity condition so is usually a very good approximation. The modes do not influence each other 
at leading order, and our calculation gives the zero-mode as 

ho{v,0 = -e^'^\ (38) 
a 

and the massive modes as 

7rm^^^ ( Yi(r7i/a) J2{m£,) — Ji{m/a) Y2{m(,) ) 

hmiv, = ' ^ , ^ e^-"" , (39) 

2J{Ji{m/a)y + {Yi{m/a)y 
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which takes the value 

hmiv, = , e^-'"" , (40) 

V'(Ji(m/a))' + (yi(m/a))' 

on the brane. 

By re-writing the equation of motion in a different way, we can see more directly that the zero-mode to leading 
order gives the same answer as four-dimcnsional cosmology. By treating one of the order e terms in the equation of 
motion as a zeroth order term, one gets 



d^h dh d^h Zdh (d^h ,d^h\ 



)+0{e^). (41) 



Setting the r.h.s. to zero for the zeroth-order calculation and performing the separation of variables in the same way 
as before gives the time evolution equation 

d'^^ 2 da 5* ,2, „ 

for the time component of the zero-mode, which is exactly the same as the equation in four dimensions. Not sur- 
prisingly, the secularity condition at first order reveals no slow-time dependence for the zero-mode. So, by slightly 
dishonest accountancy with the e terms, one can make the time dependence of the zero-mode look exactly like the 
four-dimensional mode. 

At first order in e, there is a modification in the phase due to the £-dependent term in pi|l. Evaluated on the 
brane, we have 

^-J + Ya^Ois^), (43) 

so the zeroth-order part of the solution on the brane can we written 

m/a r . /re 



^{j,{m/a)f + {Y,{m/a)y 

for the massive mode, or 



exp{±zc.,„(I-|-J)} (44) 



iexp{±^fc(r+|)} (45) 

for the zero-mode. This is in addition to the contribution of the term. 

The first-order equation in e (|37|) has an integral term on the r.h.s., which describes how a mode is sourced by 
the other modes. The matrix AAmm which is a function of s, is calculated in Appendix B. Of particular note is that 
the zero-mode cannot source any other modes, that is, A^on = 0. Massive modes, however, can decay into other 
modes. Again, this is an effect which occurs at first order in e. The effects of mode-mixing would be very small later 
in the history of the universe so one would expect all the interesting signatures coming from mode-mixing to arise 
when e is not too small, that is, when the Hubble radius is only one or two orders of magnitude larger than the AdS 
length-scale. Of course, our approximation scheme breaks down as e approaches unity. 

From studies of the evolution of perturbations during a de Sitter era in the history of the universe 0, we know 
that production of massive modes is heavily suppressed. So it is reasonable to take the zero-mode as an approximation 
for the spectrum of perturbations after the end of inflation. If the AdS length-scale, Z, is sufficiently small that our 
approximation method is valid at the end of inflation, then our analysis shows that no massive perturbations will be 
sourced during the subsequent evolution, except possibly at ©(e^). 



VII. DISCUSSION AND CONCLUSIONS 



The approximation method we have used has several advantages over the alternatives. Firstly, our coordinates are 
related to the Poincare coordinates advocated by [l5j | by a non-singular transformation and so, unlike GN coordinates 
do not have a coordinate singularity. This means that the answers we get are valid right up to the horizon, allowing 
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one meaningfully to impose initial conditions on a Cauchy surface or a boundary condition of no incoming radiation 
at the Cauchy horizon. Secondly, we have not had to specify the evolution of the scale factor, a, so the method can 
be applied quite generally in different eras of the universe and at transitions between such eras; by contrast, when 
using a GN coordinate system the time evolution has to be specified in order to find approximate solutions and, even 
then, the equations are only tractable in certain cases JJJ . Finally, it should be simple to extend the method to other 
bulk geometries where the perturbation spectrum is known. 

In particular, it is worth comparing our method to the gradient expansion method. The gradient expansion method 
works by expanding the solution in the bulk in orders of the derivatives normal to the brane. The small parameter 
is exactly the same one as used here. The two methods are quite similar, the main difference being that the method 
presented here uses a coordinate system free from coordinate horizons. Our multiple-scales method can also allow 
for the presence of very massive modes, although intuition tells us not expect such modes to be important at low 
energies. 

An important case is when the initial state for our calculation is the zero-mode solution, because this is what one 
expects from a de Sitter phase in the history of the universe 0, II] ■ Our analysis shows that the evolution of the 
solution will be the same as for perturbations in a four-dimensional FRW universe and that the massive modes are 
not sourced. This is somewhat obvious from physical intuition if these modes are viewed as Kaluza-Klein particles 
because one would expect particles to decay into less massive states, not more massive ones. The picture would be 
very different if there were sources on the brane: massive modes could be sourced by the brane-based matter. For 
simplicity, and because it is a reasonable approximation to make for the tensor part of the spectrum, we have assumed 
that there is no matter source on the brane to source the bulk metric perturbations. If one were to consider scalar 
perturbations, the brane sources could not be neglected in the same way. 

Our calculation applies strictly to the tensor component of the metric perturbations on the brane since they are 
gauge invariant, but it should be simple to apply it to the vector modes. Tensor modes are gauge- invariant, but for 
vector modes there is the added complication of gauge choice; the gauge we have used here is different from the gauges 
usually used to study four-dimensional cosmological perturbations, for example, the synchronous gauge. It should also 
be possible to extend this method to allow for a tensor matter source, such as the photon quadrupole which develops 
after last scattering. In addition, we hope, in future work, to include the effect of perturbing the brane position, so 
as to apply the method to study scalar perturbations. 
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APPENDIX A: COORDINATE CHANGE 

We consider the coordinate change between Gaussian Normal (GN) coordinates and the conformally flat bulk 
coordinates. The metrics are given respectively by 

ds2 = -iv2dr2 + dC^ + A^S^jdx'dx^ ^ — (-dT^ + dZ^ + ^^Ax'dx^) . (Al) 

We use a trick to simplify the transformation by noting that — dT is a Killing vector, so that its inner product 
with the normal vector d/dC is independent of C along the geodesies in the C direction, and hence we can write 

dT 

g^-^i^)- (A2) 

The usual formula for coordinate change gives 
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Solving these simultaneously, we find 

= N'^E'^ . (A6) 



/9Z 



On the brane, N = A = a(r) and Z = l/a so 



dZ_ 

dr 



I da , , ^ 

= = A7 



Thus, E{t) = —e/air), where the sign can be determined by considering the components of the normal vector in the 
T, Z coordinates. So, the Jacobian matrix and its inverse are 



/ dT 


dZ \ 




dr 


dT 




dT 


dZ , 




\ a< 


9C / 











)■ i 


dT 






K. 






dZ J 





Z^ , ] \ djL ^ ] \ /I /i 1 e^^^ 

+ \ dZ dz \ -7[N ^ 



(A8) 



APPENDIX B: INTEGRALS OF PRODUCTS OF BESSEL FUNCTIONS AND ORTHOGONALITY 

RELATIONS 

1. Standard results 

There is a standard formula for Bessel functions (see, for example, Ref. which tells us that 

xB^{mx)B,j(nx) dx = — {mB^^i{mx)Bi,{nx) — nBy(mx)Bu+i{nx) \ , (Bl) 

where B^ and Bjj are any two cylinder functions of order v and m ^ n. We will apply this to two specific cylinder 
function defined by 

C\r\x) = mYi{ms)J^{x) - mJi{ms)Y^{x) , (B2) 
D'^J^^x) = Y2{ms)J,{x) - J2{ms)Y,{x) , (B3) 

with J and Y being Bessel functions of the first and second kinds respectively. Cylinder functions satisfy the recurrence 
relations 0| 

zB,^i{z) + zB,+i{z)^2vB,{z) and B,^i{z) - B,+i{z) = 2'^^^ , (B4) 

dz 

from which it can be seen that 

^C("Vma;) = ^C^^^frnx) - rr? D''J^\mx) , (B5) 
OS s 

where the l.h.s. is also a cylinder function of order v. The Wronskian of Bessel functions of the first and second kinds 
is [II 

W[J„,r.] W - Jl{z)yu{z) - Yl{z)Mz) = Y,^,{z)Mz) - J,^i{z)Y,{z) = — , (B6) 
which, in conjunction with the recurrence relations, allows us to evaluate the following 

1 \ > TTS ^ OS ^ ' TTS-^ ' 

-3 ^ ' Trms'^ ' OS o \ / Trms'^ tts ' ^ ^ 

which will be of use to use later on. 

We will also make use of the asymptotic expansions for large argument 

TTX 
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the integral 

' x-P+^Bp{x)dx = -x-P+^Bp-i{x) , (B9) 



and the orthogonahty relation 



when v is a whole number. 



1 

xJy{mx)Jv{nx) — — 5{rn — n) , (BIO) 



2. Orthogonality relation 

From I)B1|) it is clear that, for ni ^ n 

X Cl"'^ (mx) Cl'''> (nx) dx = . (BU) 



because the derivatives of the C^™^ functions with respect to z vanish on the boundary z ~ s. To determine the 
coefficient of the delta function in the orthogonality relation, we note that the integral is dominated by the behaviour 
near infinity, so we can consider instead the leading order asymptotic expression 

C2 (ma;) ~ W \Yi{msf' + Ji{ms)^j sin (mx + phase angle) , (B12) 

valid when x is large. So the constant of proportionality can be deduced to be 

(nx) dx = m(Yi{miy + Ji{ml)^^ S{m — n) , (B13) 

by comparison to the equivalent asymptotic expression for eq. IjBlOjl . 

This argument breaks down in the case of the zero-mode, which must be evaluated separately. The inner product 
of the zero-mode with itself is simply 



e ° " Js e 2' 



i^s^^^, (B14) 



while (jB9|l gives us 

pOO 



^$(^)$r=0. (B15) 



3. Resonant terms 

The integrals in the first order equation of motion H25|) which need to be calculated to obtain the secularity condition 



are 



^$(0)^ (B16) 



(B17) 



To evaluate these when m = is simple, because C'q*'^ = s, and we can easily see that 
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To evaluate (|B16|I when m 7^ 0, we can differentiate the orthogonahty relation H28() w.r.t. s to get 



$(0) 



9$ 



(0) 



+ ^ '^^"^ ] - i (Mms))' + {Y,{ms)f ] 6{m ^ n) + .)<i>(°)(., .) . (B19) 



ds 



ds 



r2 

4 ""ds 



and note that, for to = ?i, the last term is ignorable compared to the delta function term. Evaluating ljB17|l when 
771 7^ 0, can be done integrating by parts to get 



^2 m 



$(0) $(0) + i 

3 771 m '2 



f-2$(0) $(0) 



(B20) 



which is just given by the orthogonality relations since, again, the second term is insignificant compared to the first. 



4. Mode-mixing terms 



The mode-mixing, that is, the effect of modes as source terms for other modes, arises from the terms involving the 
following integrals 



^ $(0) ^ 

^3 ™ ds 



Js e 



(B21) 
(B22) 



when m ^ n. These integrals are simple to calculate when either 7ti = or 77, = 0, using = s. For (|B2H) we 
calculate that 



ds 



0, 



d^^(o) 9^ 



(0) 



while IIB22I1 is evaluated to be 



dC ^(0) 5$^"^ 



e ° ds 



d£, n 



Crins) 



by applying the relation (|B9p . 

When 777, and 77 are both non-zero, we can evaluate IIB2ip using the formula (|Bip to give 



^3 ™ 



9s 4 

The recurrence formulae (IB4II allow us to write 

[4 



-Y2{ms) — 777y3(7ns) 



-J2(tos) — mJ^ims) 



d 



whence, we see that 



77^ ) Y2{ns) Y^{ns) 



(0) 



77^ ) J2{ns) Jsins) 



^3 9s 



(to + n)ms^ [m? — 772)s 



To calculate (IB22II wc first note that 



d_ 

dl 



(B23) 



(B24) 



(B25) 
(B26) 

r,(77e), (B27) 

(B28) 
(B29) 
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and similarly for ¥2- Now 



and 



2^ 



(0) 



dn 



1^1 
2 n 



d I Yi{ns) 



J2K)-n' 



dn \ n 

2 d I Yi{ns) 



d ( Jiins) 
dn \ n 



\Y2{0 



n 



dn \ n 



n dS. 2^ 

- ^^■<HY2{ns)J2{nO - J2ins)Y2{nO 
n dt, 2 V 



2 d I Ji{ns) 



dn \ n 



\Y2{nO 



so that 



and we can write 



d^^nd^ , n 



+ ^nhi(Y2{ns)J2{n^) - J2{ns)Y2{nC} 



(0) 



dn], e + 2 A e " ' ^ 



(B30) 



(B31) 



(B32) 



(B33) 



The first integral on the r.h.s. is zero when m ^ n; the second integral is 

_-2 /"OO 2 /-oo 2 /"CXD Ci 

^nhj^ ^ct\mODi-\nOd^^^J^ ^ct\mOCt\<)di~ -sj^ ^Ct\n0 (B34) 

the first term of which is zero for m 7^ n by the orthogonality relation, the second having already been calculated 
above. 
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